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Abstract. A sample road network is incrementally thinned with 34 steps 
using the Thin Road Network tool in ArcGIS® for topographic mapping 
through scales from 124,000 to 11,000,000 to evaluate procedures for use 
in The National Map of the United States. Density partitions that preserve 
local density for urban areas are combined with geometrically increasing 
scale ranges and calculations based on the Radical Law to select a reduced 
set of 12 laddered thinnings. These levels are also used to prioritize road 
labeling, so more important roads are labeled before over minor roads with 
less network connectivity. Each scale range uses four thinning levels from 
the set of 12, with the least thinned adding more local roads in the dense 
partition, the most thinned level producing roads for the sparse partition, 
and labeling prioritized with the assistance of two intermediate thinning 
levels. 
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1. Introduction 

This paper reports on progress in generalization and selective feature re- 
moval for roads in the transportation base map layer that enables compe- 
tent mapping through scales ranging from 124,000 to 11000,000. Thin- 
ning methods are applied to road features and road labels for The National 
Map of the United States, served by the United States Geological Survey 
(USGS). Roads are thinned adaptively, which removes features using their 



hierarchy within categories and network connectivity, yet preserves charac- 
teristic urban/ rural local density patterns that can be lost through simple 
category removals. For example, most U.S. roads are in the local category, 
but completely removing this category in response to urban road pattern 
congestion at smaller scales leaves some rural areas devoid of roads. The 
paper also reports progress on combining multiple levels of road thinning 
to improve the quality of road label placement (Brewer etal. 2012, in press). 



2. Background 

Road network generalization is a topic of recent interest. Recognizing the 
advantages of two network thinning approaches, Li & Zhou (2012) used a 
stroke-based approach (Thomson & Richardson 1999) in combination with 
a mesh-density approach (Chen et al. 2009) to thin road networks while 
maintaining areal patterns. Benz & Weibel (2013) extended this approach. 
Road thinning operations may be enabled or enhanced through road net- 
work or block structure (J iang & Claramunt 2004, Touya 2010, Gulgen & 
G6kg6z2011). 

Li & Choi (2002) examined road elimination through scalefor topographic 
mapping of Hong Kong at six scales from 11,000 to 1200,000. They analy- 
zed trends in road elimination related to various road characteristics, and 
found that road type was the best predictor of road elimination followed by 
length, number of lanes, number of traffic directions, width, and connectivi- 
ty at ends of road segments. The U .S. roads data used for The National Map 
does not consistently contain all of these attributes so we cannot build di- 
rectly from th ese r esu I ts. 

Gill (1993) and Phillips & Noyes (1977) noted characteristics important to 
road map reading, which are partly relevant to topographic map design. 
Phillips & Noyes found that use of straighter names (compared to strongly 
curved labels), larger sizes, and fewer size categories allowed names to be 
found faster. Gill found that people performed better on a route planning 
task if they had visual information about route continuation because they 
prefer major routes and routes with fewer turns. Participants desired in- 
formation on built-up areas so they were able to avoid towns, as well as 
avoiding hilly and winding routes. These results from older projects confir- 
med that simpler road patterns on smaller- scale displays were those that 
show and label major routes but also retained urban road densities. 

Recent research by Skiles & Howarth (2012) evaluated the advantages of 
using road signs that include simplified maps showing basic topological 
relationships among nearby towns and points of interest with the connec- 
ting roads in simplified form. Though their research is on outdoor road 



signs, it emphasizes the importance of simplified topological I y- accurate 
road patterns for spatial learning. 



3. Study Area 

Previous work used a sample area near Atlanta, Georgia (USA), to evaluate 
Esri ArcGIS® lO.lThin Road Network tools and refine their use (Brewer et 
al. 2012, Stanislawski et al. 2012). This paper reports on the next stage in 
the development of road generalization procedures intended to be useful to 
the USGS for topographic mapping through scale. It reports on results for a 
sample area surrounding Des Moines, Iowa (USA), which had a 2010 popu- 
lation of about 203,400. The study area is approximately 140 by 155 kilo- 
meters (km), with over 33,000 km of roadways (104,516 road segments; 
Figure 1). This area was chosen because it includes a variety of dense and 
sparse road patterns. 1 1 also has vari ed densiti es of hydrographi c features as 
well as corresponding flat and hilly terrain, which we are evaluating in rela- 
ted research (Buttenfield etal. 2013). 



4. Road Data Processing 
4.1. Generalization 

Two ArcGIS® generalization tools were used initially, Merge Divided 
Highways followed by Collapse Road Detail, to remove doubled and small- 
detail lines. The USGS Best Practices roads dataset used in this study inclu- 
des many dual carriageways (two parallel lines for divided highways and 
boulevards) and these were col lapsed to single lines using a 27 m tolerance 
within road classes (Figure 2a). Traffic circles and cul-de-sac ends were 
removed using a 75 m tolerance (Figure 2b), though the Col I apse Road De- 
tail tool functioned less reliably for removing small turn lanes. These gene- 
ralization tools were run on the 124,000 (24K) data and the resulting road 
lengths are listed in Table L They reduced the total road length by 882 km 
(from 33,178 to 32,296 km). 

These road details being generalized were useful for road display at large 
map scales but interfered with Thin Road Network processing at longer 
itinerary lengths. They seemed to provide many alternative routings that 
caused areas with dual carriageways to persist as high density groupings at 
smaller scales. The generalized total road length is used as the basis for sca- 
le calculations at each thinning level (discussed below; Figure 3). The un- 
generalized roads are shown at 24K-39K and for the dense partition at 40K- 
59K. The generalized roads are used for thesmaller scale ranges. Additional 



cartographic generalization such as simplification will likely be useful but 
was not applied for the results reported here. 




Figure L All roads for the Des Moines, Iowa (USA), study area at 11,350,000. 
Wide red roads are interstates and thin light brown roads are local roads. 




Figure 2. Examples of a) Merge Divided Highways and b) Collapse Road Detail 
generalization results at 130,000. Magenta roads are merged/ col lapsed to the 
lighter blue lines, enlarged in c)— ramps and intersecting roads are repositioned to 
join the single line resulting from the merge, retaining topology in the network. 



Sparse Partition 



minimum 






road length 


total road 
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length (km) 


scale 


all 
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9,723 
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8,831 


170,325 
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8,200 


197,571 
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7,747 


221,351 
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7,199 
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6,651 


300,260 
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6,159 


350,213 
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5,644 


417,040 
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5,398 


455,860 
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5,173 


496,331 
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4,772 


583,318 
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4,144 
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888,105 
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3,697 
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3,434 


1,126,133 



Dense Partition 

minimum 



road length 
(km) 


total road 
length (km) 


calculated 
scale 


all 


8,770 


24,000 
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7,176 


35,850 
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1,244 


1,193,702 


24 
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1,089 
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28 


1,048 
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1,048 


1,681,819 


30 


1,043 


1,696,029 



Table L Road lengths after generalization (using ArcGIS® Merge Divided High- 
ways and Collapse Road Detail tools) and with laddered minimum-length pro- 
cessing using theThin Road Network tool. 



4.2. Road Thinning 

We next processed the roads using the Thin Road Network tool in a ladder 
fashion, incrementing by 0.5 km minimum length up to 4 km and then byl 
km i ncrements to a 30 km minimum length, producing 34 levels (described 
here as thinning levels). A chosen minimum length parameterizes Thin 
Road Network geoprocessi ng, which makes calculations based on how 
many itineraries having lengths approximately twice this minimum para- 
meter would travel along each road segment. Table 1 shows the resulting 
total lengths of roads that remain after each successive thinning for each 
density partition. For example, the dense partition (right side of Table 1) 
starts with 8770 km, and is first thinned using a 0.5 km minimum road 
length to produce 7176 km of roads. Thinning these resulting roads with a 1 
km minimum length produces a 5999 km network, and so on until the 30 
km minimum length increment produces a road network reduced to a total 
of 1043 km. 

Ta bl e 1 al so shows the seal e at wh i ch the densi ty parti ti on i s cal cul ated to be 
suitable for display for each thinned road network. For example, when 
thinning to4 km minimum length, the sparse partition (left) is suitablefor 
display at about 132,000 scale and the dense partition is suited to 
1118,000. These scales are markedly different, and our recommended pro- 
cedures revolve around taking advantage of these differences for an objecti- 
ve approach to retaining differential densities across the map that portray 
variation in human land use. Using a single thinning level across the entire 
study area, which includes multiple land use types, homogenizes the 
appearance of the road network. As a comparison, Figure 3a shows all 
roads at 1225,000 and Figure 3b shows only the 7 km minimum road 
length thinning for the entire area (7 km for the dense partition only is 
shown at thi s seal e i n a I ater f i gure) . 

The calculated scales for each thinning level listed in Table 1 are produced 
by inverting a variant of the Radical Law that uses total length rather than 
number of features (Topfer & Pillewizer 1966; Brewer et al. 2012, in press). 
The target scales were each calculated using the total generalized road 
length at 124,000 as the source data (after merging and collapsing) and the 
resulting road length at each thinning level. The sparse and dense areas of 
the map were partitioned before beginning the road and label geoproces- 
sing (Stanislawski etal. 2012). 




a) b) 

Figure 3. Comparison maps showing a) all roads at 1225,000 and b) a more ho- 
mogeneous network resulting from roads thinned using 7 km minimum length for 
both dense and sparse partitions. 



4.3. Scale Ranges 

Our goal was to produce mapping through the range of 124,000 (24K) to 
11,000,000 (1M). To make this task more manageable for USGS, we plan- 
ned a set of nine increasing scale ranges, calculated using a geometric se- 
quence: 



• 1M and smal I er seal es 

• 630K to 999K 

• 400K to 629K 

• 250K to 399K 

• 150K to 249K 

• 100K to 149K 

• 60K to 99K 

• 40K to 59K 

• 24K to 39K 

For each scale range above, we produced one level of detail (LoD) for the 
road network that is suitable for display through the range. The particular 
beginning and end points in each range were rounded from the geometric 
calculation. This general strategy is based on our experience with these da- 
ta. Small ranges at large scales, and large ranges at small scales, with these 
approximate magnitudes were being repeatedly selected as we attempted to 
solve this problem with a variety of road datasets. This decision on range 
calculation invites later re-evaluation. 

4.4. Thinning by Scale Range 

The next step was to select a pair of thinning levels for each scale range 
from the selection of 34 levels in Table 1 For example, a 16 km thinning is 
matched with 198K for the sparse partition and a 7 km thin produced 215K 
for the dense partition, both scales falling approximately midway through 
the range 150K to 249K (Table 2). Other thinning levels also fall in this 
seal e range ( 14 to 17 km for sparse; 6 and 7 km for dense) . We went through 
a process of winnowing these choices to maximize the number of shared 
thinnings across ranges. Tables 2 and 3 list shared choices among scale 
ranges for 7 km up to 19 km thi nni ngs. Thus, rather than runni ng 34 i ncre- 
ments of the Thin Road Network tool (which took about eight hours on a 
standard desktop computer), we are able to instead process only the 12 lev- 
els listed across the columns of Table 3, offering a substantial reduction in 
processing time and improved production efficiency. 

Figure 4 shows dense and sparse partitions for part of the Des Moines 
study area (sparse partition areas are light purple). Local roads are light 
brown; state, county, and U.S. routes are light red; and interstate routes are 
dark red. The sparse area isthinned with 23 km minimum lengths (500K is 
in the400K-629K range), and the brown to red roads are present in both 
the sparse and dense thinnings. The blue roads are those selectively added 
by the 12 km thinning used for the dense partition only to provide a sense of 
greater density in the urban area without the clutter of using all roads (Fig- 
ure 3a). 



scale range 


minimum length - 
sparse (km) 


minimum length - 
dense (km) 
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29 
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150- 249 K 


16 
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100- 149 K 


12 


4 
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9 


2 


40-59 K 


7 


1 


24-39K (labels) 


0.5 


0.5 


24-39 K 


all 


all 



Table 2. Minimum lengths of thinning levels used for lines and labels in sparse 
and dense partitions within each scale range. 
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Table 3. Four thi nni ng levels were used for each scale range:three greens for dense 
partition (to left along a row) and light purple for sparse partition (right). All road 
I i nes are shown at the largest scale ( beige) 
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Figure 4. Partitions for a portion of the Des Moines study area (sparse partition 
areas are in light purple) at 1500,000. Blue roads are those added to the dense 
partition by using a shorter minimum length with the Thin Road Network tool (12 
km thinning). Reddish and brown roads arein both partitions (23 km thinning). 



4.5. Labeling Using Thinning 

We labeled local roads using a series of thinning levels (columns; Table 3) 
from the twelve- step series used across all scale ranges (rows). Dense roads 
were the I east thinned (all were shown inthe24-39K range, Figure5a) and 
label placement for the dense partition was prioritized using three thinning 
levels (label classes X 2, 3 with green tints). Sparse roads and labels were 
produced using the same thinning level (purple tint in Table 3), and these 
were one level higher from those used for the dense partition within each 
scale range. Continuing with the 150K-249K scale range as an example 
(Figure 5c), dense partition roads were thinned using a 7 km minimum 
length and three road classes at 7, 9, and 12 km were used to prioritize label 
placement. All roads in the sparse partition were labeled using the 15 km 
thinning level. 

This pattern across four levels was maintained through all scale ranges (Ta- 
bles). It is the reason the 0.5 km thinning was added, so it could be used to 
reduce the likelihood that small roads were labeled before more i mportant 
ones for the 24K to 39K scale range. This strategy is intended to put labels 
on the more i mportant roads within the local roads category for each scale 
range. 

Figure 5 shows the thinned road network results with the dense and sparse 
partitions combined on single maps through a range of example scales. 
Dangling road segments along partition boundaries were removed manual- 
ly for these figures, with that process being automated for upcoming work. 
It is important to note that these maps will eventually have other labels for 
features such as hydrography, points of interest, populated places, and 
emergency services, so smartly prioritizing road label density is as much a 
challenge as reducing road density. Likewise, the final maps will also in- 
clude more road labels as we process all names associated with a line- 
shield- labeled roads often have names in addition to number designations. 
And many road segments have multiple route numbers, not shown in these 
draft maps through scale in Figure 5. 



5. Concluding Comments 

Appropriate reduction in road network density for small-scale mapping can 
be automated if guidelines for thinning are structured in advance. These 
guidelines for successive road network thinning, however, should be gene- 
rated concomitantly with density partitioning, because using a single level 
of density or detail reduction across a dataset results in unrepresentative 
homogenization of road density, obscuring land use and morphological pat- 
terns. We have shown that a systematic reduction of density can be achie- 



ved through a relatively wide range of scales. We have also systematized 
label prioritization through scale. Future work will involve the application 
of these guidelines to other datasets with different road network patterns. 
Wewill also evaluate label hierarchies within the road classes. 
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Figure 5. Example scale series: a) 130,000; b) 175,000; c) 1225,000; d) 
1500,000; and e) 11000,000. These maps are thinned and labeled according to 
Tables 2 and 3. 




6a) All roads; labelsthinned using0.5to 2 km minimum lengths. 




6b) Roads and labelsthinned using 2 to 9 km minimum lengths. 




5c) Roads and labelsthinned using 7to 16 km minimum lengths. 




5d) Roads and labels thinned using 12 to 23 km minimum lengths. 




5e) Roads and labels thinned using 19 to 29 km minimum lengths. 



